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Abstract. A general framework for calculating lunar and solar torques on the oceanic tides is 
developed in terms of harmonic constituents. Axial torques find their associated angular momentum 
and earth-rotation variations are deduced from recent satellite-altimeter and satellite-tracking tide 
solutions. Torques on the prograde components of the tide produce the familiar secular braking of 
the rotation rate. The estimated secular acceleration is approximately -1300'ycentury 2 (less 4% after 
including atmospheric tides); the implied rate of change in the length of day is 2.28 milliseconds/century. 
Torques on the retrograde components of the tide produce periodic rotation variations at twice the tidal 
frequency. Interaction torques, e.g. solar torques on lunar tides, generate a large suite of rotation-rate 
variations at sums and differences of the original tidal frequencies. These are estimated for periods 
from 18.6 years to quarter-diurnal. At subdaily periods the angular momentum variations are 5 to 6 
orders of magnitude smaller than the variations caused by ocean tidal currents. 
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1. Introduction 

Brosche and Seiler (1996) recently called attention to the interesting role that direct lunar and solar 
torques on the ocean tide play in the earth’s short-period angular momentum balance (“short-period” 
here meaning daily and subdaily). Brosche and Seiler noted that such torques are a potential source 
of angular momentum in the earth-ocean system and that this source had been neglected in previous 
earth-rotation studies. The purpose of the present paper is to reexamine, clarify, and extend these 
ideas. We limit the discussion to the earth’s spin rate and ignore the additional complications of wobble 
and nutation. It suffices therefore to examine only axial torques. We also limit the discussion to diurnal 
and semidiurnal tides of the second degree in the tidal potential, i.e. to the “major” short-period tides. 

A qualitative understanding of the consequences of tidal torques can be obtained from the diagram 
in Figure 1, which is drawn for the principal semidiurnal tide M 2 . Two spherical harmonic components 
of the ocean tide — the only two that induce nonzero torques — are displayed: prograde and retrograde 
components of degree 2, order 2. The prograde component is the classical tidal “bulge” that propagates 
westward under the moon; the retrograde component is a similar bulge, generally smaller, that 
propagates eastward. (“Prograde” here implies moving in the same direction as the tide-generating 
body; this prograde /retrograde terminology follows standard tidal usage, but it is opposite that used in 
polar-motion studies.) The retrograde bulge owes its existence, of course, to the highly nonequilibrium 
form of the oceanic tides; there is no comparable bulge in the body tide. The tidal force acting upon 
the prograde bulge produces the familiar secular braking torque (Munk and MacDonald, 1960). The 
tidal force acting upon the retrograde bulge produces a periodic torque at twice the tidal frequency. 
This latter phenomenon is quite similar to that of libration (Chao et al, 1991; Wiinsch, 1991) in which 
the triaxial ellipsoid of the solid earth experiences a periodic torque at exactly the tidal frequency as 
it rotates beneath the moon. But because the ocean’s retrograde bulge travels eastward at twice the 
earth’s rotation (relative to the moon), the torques on it are twice as rapid. 

In the same manner, there are interaction torques between tidal constituents — for example, solar 
tidal torques on lunar tidal bulges. These torques have periodicities at the sums and differences of the 
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original constituent frequencies, the summed frequencies corresponding to torques on retrograde bulges 
and the differenced frequencies corresponding to torques on prograde bulges. These new frequencies are 
actually those of the familiar compound tides that occur in many shallow seas when nonlinear dynamical 
terms (e.g., quadratic bottom friction) cause interactions between tides. (The two phenomena are, 
of course, unrelated.) We thus find a whole new set of frequencies to explore in the earth’s angular 
momentum budget, previous studies having been confined to only the tidal frequencies. Intriguing as 
they are, however, these new terms turn out to be quite insignificant. 

In the following calculations, earth anelastidty is ignored. This allows us, among other things, to 
neglect the torques on the lagged body tide and to concentrate on the oceanic problem. Given the 
most recent estimates of the body-tide lag (e.g. Ray et al., 1996), neglecting anelastidty may induce 
errors of perhaps a few percent at most. There are, nonetheless, still nonzero torques on the solid-earth 
tide, since there is a fairly large radial-displacement load tide caused by the oceanic tide, which acts to 
reduce the overall torque by some 30%. 

Finally, it is important to realize that the torques under discussion are on the tidal height 
fluctuations alone, not the entire ocean. Torques on the latter, integrated vertically from mean sea level 
to the ocean depths, are included in the solid-earth libration, which is dependent on the inertia tensor 
of the whole earth, including the ocean, and is determined from the gravitational Stokes coeffidents of 
the whole earth. This point is indeed stressed by Brosche and Seiler (1996), but it is undear from their 
paper whether they followed it, since they found all torques of the same tidal spedes to have precisely 
the same phase (see JT' in their Table 1). A constant phase does occur for libration, but it should not 
be expected for torques on the ocean tides. 

2. Torque Estimates 

Much like the precession/nutation problem (Moritz and Mueller, 1987, p. 52), there are two 
ways to compute the required torques. One may compute the direct tidal torque on the ocean tide 
(and its load deformation), or one may compute the equal and opposite torque of the ocean tide’s 
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gravitational attraction (including its load deformation) on the moon (or sun). The advantage of the 
former method is that one may use directly the readily available harmonic developments of the tidal 
potential (e.g., Cartwright and Tayler, 1971) and thereby avoid dealing with the complications of the 
lunar and solar ephemerides; the disadvantage is that the torque on the load deformation, requiring 
volume integration, is difficult to calculate directly. We nonetheless prefer this method, because by 
using the Cartwright-Tayler expansion we find it easier to sort out the frequencies and phases of torque 
components and their underlying ocean constituents. In addition, the complicated torque on the earth’s 
load deformation may be handled trivially by exploiting the symmetry of the two torque methods, as 
will be seen. 

Note that for the problem at hand, torques internal to the ocean-earth system may be ignored. 
They would include, for example, the mechanical force of the body tide against the ocean tide and the 
Newtonian attraction of the ocean tide on the body tide. In the present analysis, such internal torques 
serve only to couple the ocean and solid earth into a mutually rotating solid body. 

From the introductory discussion, it is evidently advantageous to express the ocean tidal height, 
for any given harmonic constituent of frequency <Ti, in the form (Lambeck, 1980, eq. 6.2.1 — allowing for 
a misprint in 6.2.1c) 

C = X) X) Y! D nm COs(<7if ± TTUj> - V^ m ) P n m (cOS0), (1) 

n m + 

comprising prograde (+) and retrograde (-) waves of amplitude D± m and phase lag where 

(0, (j>) are spherical polar coordinates and an associated Legendre function. Lambeck’s notation 

D ± cos(a ± 0) denotes D+ cos(a + 0) + D~ cos(a - 0). This C is the ocean’s tidal height fluctuation 
relative to the seabed. The D± m and tp^ m parameters are known either from numerical hydrodynamic 
models of the global tide or from space-geodetic measurements (e.g. Schrama and Ray, 1994). 

Following Cartwright and Tayler (1971), the astronomical tidal potential at the earth’s surface for 
a constituent of frequency (72 is written 


U = gH ^(5/2471-) (cos 9) si n(a 2 t + 4>) 


(2a) 
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for diurnal tides and 

U = gHy/(5/96n) P 2 2 (cos 6) cos(<r 2 * + 2<f>) (2b) 

for semidiurnal tides, where g is the gravitational acceleration and the amplitudes H (in length units) 
are as tabulated by Cartwright and Edden (1973). Note that frequency <j\ refers to the tidal elevation 
(1), while cr 2 refers to the potential (2); each may be either diurnal or semidiurnal, and ci may or may 
not equal <72. 

The potential U acting on the ocean tide £ produces an axial torque of general form 

r;= / p(xx vu) z dv, 

Jv 

where p is the density of seawater and x is a vector from the geocenter to the differential mass element 
pdV. This torque can be simplified to a surface integral by the thin-shell approximation: 

r 1 = J pr sin 9 (VU)# dS = p J C (dU/d<f>)dS, 

with dS a differential of area. The corresponding torque of U acting on the solid-earth load tide cannot 
be similarly simplified, but it can be trivially accounted for, as noted above, by considering the torque 
of the ocean tide’s gravitational potential acting on the moon, which is equal to — Y° z . For that torque, 
accounting for the load deformation requires only an additional factor of (1 4* fcj), where k 2 is a load 
Love number (Munk and MacDonald, 1960). The equivalence of the two torques implies that T z must 
also include the torque on the deformation via the identical factor. Hence the combined torque is 
written 

r z = (1 + k‘ 2 )p ! C (dU/d4>) dS = { 1 + k ' 2 )pal f C (8U/dj) sin 6d0d<j>. (3) 

Js Jo 

With the previous expressions for £ and {/, integration is straightforward. Owing to orthogonality 
of the spherical harmonics, only terms in ( with (n, m) = (2,1) for diurnal tides and (2,2) for semidiurnal 
tides are effective. (An identical simplification occurs for work and dissipation integrals; see, for 
example, Lambeck [1980] and Platzman [1984].) Inserting (1) and (2) into (3) and expanding terms, 
one finds intermediate integrals of two types (note that m, l = 1 or 2 and that S m £ is Kronecker’s 
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(5-function): 


and 


r2ir 

I sin(ayt 4 m0) sin (<72* 4 t<f>) d<j> — it 5 ml cos(ai* — atf) 
Jo 


/•2tt 

/ sin(<Jit + m<j>) sin(<72* — f<^) = — 7r^m/ cos(<7i* 4 cr2*)> 

Jo 

and similar terms involving cosines and sine-cosine products, each of them yielding either sums or 
differences of tidal frequencies. The total resulting torques are 


r z = -y/(6n/5) (1 4 k' 2 ) pGMH^^D^i cos(ayt 402* ±^21) 


(4a) 


for diurnal tides and 


T z — -^(96^/5) (1 4 A^) pGMHy^Pf 2 sin(<7i* 4 a^t ± tp^) (4b) 

+ 

for semidiurnals. The earth’s gravitational constant GM has been used in place of gd When a tidal 
constituent’s potential is acting on its own tidal bulge, then ay = a 2 (h a) and the torque integrals 
simplify to 

T z = — v/( 67 t/ 5) (1 4 k^pGMH cos ^21 4 D ^ cos(2a- — ^21)} 
for diurnal tides and 


r, = -v/(96tt/5) (1 + kfipGMH sini/>+, + D~ 2 sm(2a - V> 2 ~ 2 )} 

for semidiurnals. In each expression the first term on the right is the secular braking torque. 

The implications of T z for the earth’s angular momentum budget are considered below for three 
frequency regimes: secular (<7i — a 2 = 0), long period (ay — a 2 / 0), and short period (ay 4^). This 
requires estimates of the ocean-tide coefficients £>nm>V'nm» to which we now briefly turn. 


3. Oceanic tide coefficients 

For the oceanic tide coefficients we adopt (see Table 1) three sets of estimates determined from 
satellite geodetic measurements, one from Geosat altimetry, one from Topex/Poseidon (T/P) altimetry, 
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and one from multiple-satellite tracking. In the case of altimetry, the coefficients are 

determined by numerical quadrature of deduced global oceanic cotidal charts. In the case of satellite 
tracking, they are determined directly from the tidally induced orbited perturbations (Lambeck, 1980), 
often as a part of a simultaneous, large-scale inversion for the earth’s gravitational Stokes coefficients 
(e.g., Christodoulidis et al., 1987). The T/P and satellite-tracking solutions are more accurate than 
any available numerical hydrodynamic model of the tides; see Shum et al. (1997) and Figure 1 of Ray 
et al. (1996). 

Because altimetry is a geometical measurement while tracking is a gravitational measurement, 
there are some important differences. For example, they differ significantly in estimates of the S 2 tide, 
since the tracking is also sensitive to the S 2 atmospheric tide (e.g., Cartwright and Ray, 1991); a smaller 
effect may also occur in Ki owing to seasonal variations in the S\ atmospheric tide. Both altimeter and 
tracking solutions also differ very slightly owing to the different ways that they are affected by earth 
anelasticity (Ray et al., 1996). These small anelastic effects may be ignored for present purposes. 

The Geosat satellite altimeter determination is taken from Cartwright and Ray (1991). These 
tidal coefficients are less accurate than those from the later T/P mission (e.g., Shum et al., 1997), and 
they should now be considered obsolete. They are used here to compare deduced secular acceleration 
estimates to previously published estimates from the same model. 

The T/P tidal coefficients are updates to those published by Schrama and Ray (1994). In addition 
to incorporating more altimeter observations and improving horizontal resolution, this solution uses 
a half-dozen supplemental hydrodynamic tidal models in various marginal and polar seas to ensure 
complete coverage of the global ocean. The hydrodynamic model of Le Provost et al. (1994) is used for 
all latitudes above the 66° turning latitude of the T/P satellite. The M 2 estimates of £> 22,^22 were 
previously published in Ray et al. (1996). 

All satellite altimeter estimates of ocean-tide coefficients have been corrected for the geometrical 
effects of the body tide and the load tide. The latter is based on the iterative method outlined in 
Cartwright and Ray (1991; appendix A) and used loading Love numbers h! n from Farrell (1973). For 
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T/P, the body tide was computed with Love number /12 = 0.609 for all tides except Ki and its nodal 
sidelines for which = 0.52. For Geosat, /12 = 0.619 was used in the original mission data products, 
but our M 2 and S 2 coefficients were later adjusted to conform with h% = 0.609 (Cartwright and Ray, 
1991). 

The third set of estimates is taken from EGM96S (Lemoine et al., 1998), a satellite-only solution 
for the earth’s gravitational field based on optical, doppler, laser, and GPS tracking of 40 artificial 
satellites. Only tides that cause significant long-period perturbations in satellite orbits were included in 
this solution, so they are limited to prograde terms only (Lambeck, 1980). 

The EGM96S model provides well-calibrated estimates of tidal coefficient uncertainties, for which 
the reader may refer to Lemoine et al. (1998). For M 2 the EGM96S standard error for Djj sin 1^22 13 
0.024 cm; other semidiurnal tides are comparable, while diurnal tides tend to be 3 to 4 times larger. 
Unfortunately, owing to the manner in which they were developed, standard errors for the two altimeter 
solutions are not available. However, M 2 errors have been estimated for two other T/P solutions — by 
Ray et al. (1994) and by Egbert et al. (1994)— and these should be comparable to the Schrama and Ray 
(1994) solution; for M 2 these uncertainties are 0.036 and 0.026 cm, respectively, which are comparable to 
EGM96S. Work is in progress (G. D. Egbert, personal comm., 1998) to establish a more comprehensive 
error model for the altimeter tide solutions. For our purposes, the given rough figures are adequate. 

4. Secular braking of the earth’s rotation 

In this section, we examine only those terms in T z having no explicit time dependence, taking 
ci = <r 2 . These terms, dependent only on the prograde oceanic tide coefficients, completely determine 
the present-day secular braking of the earth’s rotation rate. That secular acceleration is given by 

Q = T Z /C 

where T z here represents the appropriate relevant terms in (4a, b). Here C — 8.036 x 10 37 kgm 2 is the 
polar moment of inertia of the whole earth. 

Estimates of the secular terms from the three tide models are listed in Table 2, along with 
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standard error estimates from EGM96S, The potential amplitudes H and the (slightly) frequency- 
dependent factor (1 4* k 2 ) are also tabulated; the former are from Cartwright and Edden (1973), the 
latter from Wahr (1980). Other required constants are: 
p = 1035 kgm“ 3 
GM = 3.986 x 10 14 m 3 s” 2 . 

As is well known, the dominant tide in SI, responsible for about 70% of the total acceleration, is 
M 2 . The next largest tide is S 2 , but its oceanic effect is partly cancelled by the S 2 atmospheric tide, 
which tends to accelerate the earth’s rotation rate by roughly +55" cy~ 3 (based on spherical harmonic 
coefficients listed in Haurwitz and Cowley, 1973). The totals listed in Table 2 are misleadingly close, a 
coincidence of favorable cancellations. 

The largest tide neglected in Table 2 is probably ^ 2 * It can be estimated from the semidiurnal 
admittances; for the EGM96S model, we get D 22 = 0.121 cm, tp+ 2 = 117°, which implies fi = —1.4" cy~ 2 . 
A handful of other minor tides would add a few more arcseconds. 

In passing, we note that the long-period tides, with zonal potential 

U = gHy/(5/4n) P 2 (cos#) cos crt y (2c) 

add nothing to the earth’s secular rotational braking. With perfect symmetry in the <j> direction (dU/d<f> 
= 0), its axial torque is identically zero. This has not been sufficiently dear in some of the recent 
literature. 

Allowing for an atmospheric contribution of +55"cy“ 2 , we have a total secular acceleration of 
approximately -125G" cy“ 2 . The corresponding rate of change in the length of day is 

A = = 2.28ms/century. 

(Strictly, the units are ms/day/cy, but the day - " 1 is generally assumed.) The observed secular 
acceleration (Stephenson, 1997) includes, of course, nontidal effects as well, the most important being 
an additional +300" cy“ 2 usually attributed to melting of polar ice and/or the earth’s viscous rebound 
(Yoder et al., 1983). 
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Finally, it is of interest to compare the estimates for the Cartwright-Ray model to previously 
published estimates for that model in Ray (1994). The previous values were computed, following 
Lambeck (1977), as follows: the expression for the gravitational attraction of the ocean tide on the 
moon is converted from spherical coordinates into lunar Keplerian elements and then inserted into the 
Lagrange planetary equations (Brouwer and Clemence, 1961), yielding secular rates for the moon’s 
mean motion n, eccentricity e, and inclination I. Invoking conservation of angular momentum then 
gives the earth’s Cl. The two approaches to f l agree to within about 1% for M 2 and within 1" cy“ 2 
(the resolution of Table 2) for ail other tides. These can likely be further reduced by closer attention 
to consistency in geophysical and astronomical constants and models; this will become more critical as 
observational accuracies continue to improve. 

5. Long-period variations 

The same prograde oceanic terms responsible for the secular acceler ation also induce long-period 
oscillations in the earth’s rotation rate when <j\ ^ 02 * We consider two cases: (1) a constituent acting 
upon its nodal sideline, which generates 18.6-year variations and (2) the case when G\ y a<x correspond to 
the largest tides, M 2 and S 2 , which generates variations at the frequency of the near-fortnightly MSf 
tide. 

5.1 Variations at 18.6 years 

All lunar tidal constituents undergo an 18.6-y nodal modulation as the lunar orbital plane precesses 
around the earth. As a consequence, any tide’s contribution to the secular braking torque must undergo 
a corresponding periodic variation. Oi has one of the largest nodal modulations, so we concentrate 
initially on it. The Oi modulation in H will be augmented by other lunar diurnal tides, all of which 
are largest when the lunar orbit is maximally inclined to the earth’s equator. This occurs when the 
longitude of the moon’s ascending node is aligned with the equinox (i.e., when the longitude of the node 
is zero). These diurnal-tide modulations are offset, however, by the semidiurnal tides, which are largest 
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when the lunar orbit is closest to the equator (K 2 is an exception). As will be seen, the cancellation 
between diurnals and semidiurnals turns out to be nearly complete. 

Consider Eq. (4a) when a x is the frequency of the primary Oi line, 545.555 in Doodson’s 
nomenclature, and <72 is the frequency of the primary nodal line, 545.545. Then H = 4.945 cm 
(Cartwright and Edden, 1973). Adopting £> 21,^21 ^* om the EGM96S solution yields 

-V(6tt/5) (1 + k' 2 )pGMHDl \ cos {at + ^ 1 ) = -7.45 x 10 14 cos {at + 46°) N m 

for the relevant component of the torque IV Here, at = a x t - arf = N't = N' where AT' is the negative 
of the longitude of the lunar node (this being the argument associated with the fifth Doodson number) 
and N' the corresponding angular frequency, 1.070 x 10“ 8 s” 1 = 1 cycle/18.6 y. 

Consider now the complementary torque obtained by switching a x and a%- In principle, evaluation 
of this torque requires knowledge of the nodal line’s D% x and rp 21 . But because the deep-ocean tided 
admittance is known to be smooth (Munk and Cartwright, 1966), these parameters are equivalent to 
the primary line’s parameters, but with amplitude scaled by the ratio in H\ the known exceptions 
to this in tide-gauge data occur only in shallow, nearly resonant seas (e.g., Ku et al., 1985) and are 
unlikely to affect a global coefficient of spherical harmonic degree 2. Hence this second torque is 

-V0W5) (1 + k 2 )pGMHD 2l cos(-at + ip 2l ) = —7.45 x 10 14 cos(-c7* + 46°) Nm 
and the sum of the two torques is 

— 2v^(67t/5) (1 + k^pGMHD^ cos i) 21 cos at = -1.05 x 10 15 cosAf' Nm 

At the 18.6-y period, we again assume that C represents the moment of inertia of the whole earth. 
The amplitude of the acceleration is thus -1.31 x 10 -23 s~ 2 . The corresponding increment in the 
earth’s rotation, expressed in terms of variations in Universal Time, is 

AUT = - J AA dt = (1/n) JJ ildt 2 

= (jV 2 n ) _1 x 1.31 x 10 -23 cos N 1 


1.57 cos N* ms. 
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In similar fashion, the major lunar tides generate the following terms: 

Qx : 0.07 cos A/' ms 

Oi: 1.57 cos N* ms 

Ki: 1.90 cos AT' ms 

N 2 : -0.17 cos A 7 ms 

M 2 : -4.00 cos A/ 7 ms 

K 2 : 0.36 cos A 7 ms 

The signs of these terms are consistent with the nodal modulations of these tides, which act to reduce 
the amplitudes of M 2 and N 2 when N* = 0 and to increase the other four (Doodson and Warburg, 1941, 
Table 7.3). Of the semidiurnals, K 2 is anomalous because it is a M declinationaT tide: it arises from the 
modulation of the principal tide M 2 caused by the twice-monthly excursions of the moon away from the 
earth’s equator. This modulation, and hence the amplitude of K 2 , therefore increases, not decreases, 
with the lunar declination. 

The large cancellations between diurnals and semidiurnals leave a residual of only 0,27 ms. This 
variation in UT is very small compared to the dominant tidal variation at the nodal frequency, which 
is due to the body tide’s direct perturbation to the moment C . According to Yoder et al. (1981), the 
body tide generates oscillations in UT of amplitude 172.05 ms. 

There is no dynamical reason for supposing that the above cancellations might in fact be perfect, 
since they depend on the tidal behavior of the oceans. One might imagine, for example, an ocean that 
somehow suppresses one tidal species, causing far less cancellation than now exists. 

5.2 Variations at MSf 

In a completely analogous manner, tidal interactions between constituents generate a large suite 
of periodic oscillations in the rotation rate with periods correspondingly shorter than 18.6 y. One of 
the largest is presumably the torque of M 2 on S 2 , and vice-versa, which generates oscillations at the 
frequency of the near-fortnightly tide MSf (period 14.77 days). Because the potential amplitude of MSf 
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is relatively weak (compared to, say, the nearby tide Mf), it exists in the ocean primarily as a nonlinear 
compound tide confined to shallow water. Its direct effect on UT through the body-tide variation in 
C is correspondingly weak. The ocean-torque mechanism is therefore more important for MSf than 
it is for most other long-period tides, and in fact it might conceivably explain the anomalously large 
observed variation in UT1 at this frequency (Chao et al., 1995). The question has been examined 
recently by Cartwright (1997), who finds the effect too small to explain the anomalous observations. 
The following calculation repeats Cartwright’s in the context of the above formalism. 

Unlike the nodal variation, the lag ^ f° r ^2 and S 2 are different, and the two torques cannot be 
trivially combined as in the nodal modulation case. The variation in UT1 is (with obvious subscripts) 

AUT = {a 2 nC)~ V(96tt/5) (1 + k^pGM { H S D M sin (-<rt + rp M ) + H M D s sm(at + ^s)} 

Here a is the frequency of MSf, and at = 2D where D is the mean lunar elongation. At this frequency, 
the mantle is presumably uncoupled from the core and C should be decreased by 12% (Yoder et al., 
1981), After some algebraic manipulation, with M 2 and S 2 coefficients from the Schrama-Ray model, 
this becomes 

AUT = 0.23 cos (2D + 3°) /is. 

The excess length of day is 

A A = 0,098 sin(2D 4* 3°) /is, 

very close to that estimated by Cartwright (1997). According to Yoder et al. (1981), the main 
body-tide effect at MSf causes a variation in UT of amplitude 78.1 /is. Hence, as Cartwright noted, the 
ocean-torque interaction is of no great significance. 

6. Short-period variations 

We come finally to the subject of Brosche and Seiler (1996) who suggested that torques on the 
ocean tide may be important at daily and sub-daily periods. The torques involving Oi, Ki, and M 2 
are likely to be of comparable orders of magnitude. It suffices here to examine the torque of M 2 on its 
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own retrograde oceanic bulge, which produces variations at twice the M2 frequency. (Note that other 
combinations of tides will produce variations at precisely the M 2 frequency, for example K x +Oi, and 
M4 - M 2 . These torques are likely comparable or smaller than the M 2 +M 2 torque.) 

The angular momentum variation associated with the torque of M 2 on its retrograde bulge is (with 
or now the frequency of M 2 ) 

j = J r z dt 

= (2a) -1 /(967T/5) (1 + k' 2 )pGMHD 22 cos(2at - if> 22 ) 

= 3 x 10 19 cos(2 at - 26°) kg m 2 s -1 , 

based on the ^ Schrama-Ray estimates. This is 5 to 6 orders of magnitude smaller than the 

angular momentum exchanges between the ocean find mantle that are due to tidal currents, according 
to estimates published by Seiler (1991) and Chao et al. (1996). We conclude that for the short-period 
angular momentum balance, direct torques on the ocean tide are of no significance. 

This conclusion should be no surprise if we compare the oceanic torque mechanism with the 
closely related solid-earth libration. Libration is caused by a 70-meter difference between the semimajor 
and semiminor equatorial axes (Torge, 1980), and the torque on this solid-earth “bulge” causes UT1 
variations of order 1 /xs (Chao et al., 1991; Wunsch, 1991), which is an order of magnitude smaller than 
the variations induced by oceanic tidal currents. In contrast, the ocean tide is a bulge of only a few cm 
of water. The corresponding UT1 variations must therefore be 4 to 5 orders of magnitude smaller than 
the libration effect. 

7. Summary 

A general framework for computing lunax and solar tidal torques on the oceanic tides allows both 
the secular acceleration and a large suite of high-frequency periodic accelerations in the earth’s rotation 
to be estimated. Estimates of the torques and accelerations require only the known amplitudes in 
a harmonic development of the tidal potential (e.g., Cartwright and Tayler, 1971) and the degree-2 
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spherical harmonic coefficients of the oceanic tided elevations. The latter are available from numerical 
hydrodynamic models, from satellite altimeter analyses (e.g., Schrama and Ray, 1994), or from satellite 
orbit-perturbation analyses (e.g., Christodoulidis et al., 1987). For the secular acceleration we find 
(Table 2) a total oceanic contribution of approximately -1300"cy“ 2 . All the short-period variations 
turn out to be small. The torques on the retrograde bulges, which generate rotation variations at twice 
the original tidal frequencies, are especially small and can be ignored in studies of the earth’s angular 
momentum balance. 

Acknowledgments. We thank Drs. Shailen Desai and David Cartwright for comments on the original 
manuscript. 
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Figure 1. Geometry of the prograde and retrograde oceanic tidal bulges and the moon, seen from above 
the earth’s north pole, in a frame rotating with the earth. The prograde bulge is shown lagging the 
moon by an angle t/^/2 ~ 65° according to the estimate of Schrama and Ray (1994). The rotation rate 
is H — 7.2921 x 10~ 5 s -1 ; the moon’s mean motion is n = 2.6653 x 10~ 6 s" 1 . 



21 


Table 1. Estimates of Tidal Coefficients 


Tide £>2m, cm V’s'm- de S D 2m> cm VW de S 


Geo sat model - Cartwright & Ray (1991) 


Qi 

0.49 

39.3 

0.27 

201.5 

Ox 

2.33 

42.7 

1.12 

210.2 

Px 

0.86 

47.1 

0.47 

225.5 

Ki 

2.56 

47.2 

1.45 

225.9 

n 2 

0.70 

127.4 

0.11 

351.7 

m 2 

3.45 

131.6 

0.59 

22.5 

s 2 

1.16 

126.2 

0.31 

94.0 

k 2 

0.31 

125.8 

0.08 

99.5 


T/P model - 

Schrama & Ray 

(1994, updated) 

Oi 

2.592 

46.09 

1.160 

209.58 

Kx 

3.001 

42.43 

1.528 

224.92 

N;i 

0.702 

119.51 

0.117 

10.74 

m 2 

3.231 

129.38 

0.585 

26.07 

s 2 

1.241 

131.91 

0.217 

80.04 


EGM96S model - Lemoine et al. (1998) 
Q! 0.587 37.93 


Ot 2.732 45.09 

P L 0.991 43.59 

Ki 2.828 39.45 

N 2 


0.639 


114.97 



T^ible 1. (continued) 


Tide 

D t m > cm 

V’+m. deg D 2m , cm ip 2m , deg 

m 2 

3.266 

128.21 

s 2 

0.785 

145.90 

k 2 

0.273 

122.03 


m = 1 for diurnal tides, 2 for semidiurnal tides. 
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Table 2. Estimates of Secular Acceleration Cl by Oceanic Tides 


Tide 

H, m 

(1+ k' 2 ) 

C-R 

S-R 

EGM96 

Qi 

0.0502 

0.689 

-3 

-3 

— 3 ± 1 

0! 

0.2622 

0.689 

-63 

-66 

-71 ± 3 

Pi 

0.1220 

0.700 

-10 

-13 

-13 ± 2 

Ki 

0.3687 

0.730 

-96 

-122 

-120 ± 7 

n 2 

0.1210 

0.692 

-38 

-42 

-40 ± 2 

m 2 

0.6319 

0.692 

-924 

-894 

-919 ± 9 

s 2 

0.2940 

0.692 

-156 

-154 

-73 ± 5 

k 2 

0.0799 

0.692 

-11 

-12 

-10 ± 1 

Total 



-1305 

-1306 

-1304 ± 12 


Units of Cl are arcseconds per century 2 , 1" cy“ 2 = 4.868 x 10“ 26 s~ 2 . 
C-R = Cartwright and Ray (1991); S-R = Schrama and Ray (1994); 
EGM96 = Lemoine et al. (1998). S-R estimates of Qi, Pi, Ka are inferred 
from the admittance defined by its major tides. The EGM96 estimate of 
S 2 includes effects of the S 2 air tide, for which +55 has been added to its 


total. 





